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Abstract 

We consider the interaction of a heavy quark-antiquark pair moving in A/" = 4 SYM 
plasma in the presence of non- vanishing chemical potentials. Of particular importance is 
the maximal length beyond which the interaction is practically turned off. We propose 
a simple phenomenological law that takes into account the velocity dependence of this 
screening length beyond the leading order and in addition its dependence on the R-charge. 
Our proposal is based on studies using rotating D3-branes. 



* Currently serving in the Hellenic Armed Forces. 



Contents 

1 Introduction and summary 1 

2 Moving qq pairs in the J\f = 4 plasma 5 

3 Screening in the zero— R-chcirge case 9 

3.1 Motion perpendicular to the qq-pair axis 10 

3.2 Motion parallel to the qq-p&ir axis 12 

4 Screening in the case of nonzero R-charges 13 

4.1 Non-extremal rotating branes 13 

4.1.1 Two equal nonzero angular momenta 14 

4.1.2 One nonzero angular momentum 15 

4.2 Velocity and R-charge dependence of the maximal length 15 

4.2.1 Two equal nonzero angular momenta 16 

4.2.2 One nonzero angular momentum 18 

5 Plasma rest frame and averages 19 



1 Introduction and summary 

Since the early days of the AdS/CFT correspondence [1], there has been considerable 
interest in applying string-theoretical techniques for the computation of the interaction 
potential of an external heavy qq pair, as given by the expectation value of a Wilson 
loop, first demonstrated in [2] for the conformal case. In extensions of this work to non- 
conformal cases, it was often found that there exists a maximal length beyond which the 
interaction is practically turned off, the archetypal example being finite-temperature jV = 
4 SYM theory, described by non-extremal D3-branes [3]. This phenomenon of complete 
screening persists [4] in the J\f — A theory with chemical potentials, modelled using 
rotating D3-brane solutions, as well as in the Coulomb branch of the theory, modelled 
by continuous distributions of D3-branes, even though in that latter case half of the 
maximal supersymmetries are preserved. Although such results would be hard to obtain 
using conventional field theoretical methods, their relevance for real QCD was obscured 
both by the major qualitative differences between QCD and SYM theories and by the 
fact that the quark-gluon plasma (QGP) used for the experimental study of thermal 
QCD was until recently thought to be a weakly-coupled collection of quasiparticles [5] . 

This situation changed drastically after the realization that the QGP is a state of 
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matter best described as a strongly-coupled liquid, but still exhibiting partial deconfine- 
ment of color charges (see e.g. [6] for reviews). These features motivate the study of 
various quantities of physical interest within the gravity/gauge theory correspondence. 
The best-studied examples encouraging this line of thought are the calculation of the 
shear viscosity in A/" = 4 SYM [7] and the various calculations of the rate of energy loss 
of partons moving through the plasma [8]- [16], both of which lead to results that are in 
good qualitative agreement with experiment. 

This renewed interest shed light into some, often overlooked, aspects of Wilson-loop 
calculations of heavy-quark potentials. One such aspect is that, in real experimental 
situations, neutral qq pairs are produced with some velocity relative to the plasma, a fact 
that affects the interaction potential of the pair. Such considerations apply, for example, 
to the effect of charmonium suppression, the dissociation of J/%1) mesons due to color 
screening by the deconfining medium [17, 18], an effect that is expected to be enhanced 
when the J/%1) has a nonzero velocity with respect to the plasma. Although dynamical 
processes of this kind would be hard to study using conventional methods, they can be 
addressed rather easily within the gauge/gravity correspondence. In particular, the non- 
extremal brane metrics serving as gravity duals of finite-temperature gauge theories are, 
unlike their extremal counterparts, sensitive to boosts in the directions along the brane 
where the gauge theory resides. Then, the computation of the potential interaction of a 
qq pair moving with respect to the thermal plasma is simply accomplished by applying 
the standard AdS/CFT calculation of Wilson loops to such boosted backgrounds.^ The 
first examples of finite-temperature computations of this kind were done for a confining 
non-supersymmetric theory in [20] and for jV = 4 SYM at vanishing R-charge in [21, 
22] (see also [23, 24] for further works). The results of [21], in particular, suggest a 
universal law for the behavior of the maximum screening length was suggested, essentially 
stating that it depends mainly as (1 — v^Y^^ on the pair's velocity, perfectly valid for 
ultrahigh velocities. Given the importance of such a law for building up phenomenological 
models for certain effects occurring in hot plasmas, it is interesting to ask how this law is 
affected when chemical potentials are turned on and what are the details of the velocity 
dependence. In this paper, we examine precisely these two issues, within the AdS/CFT 
correspondence, using rotating D3-brane solutions. In principle, the maximal screening 
length can be a function of the plasma parameters, namely the temperature and the 

^To appreciate this, note, for instance, that a computation of the potential for a heavy qq pair moving 
in a hot plasma has been done only for the case of an abelian gauge theory [19]. 
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R-charge density and in addition it depends on the velocity vector of the qq pair. 

Our studies suggest a law that can be used in phenomenological models, namely that 
the maximal length has the following form in the rest frame of the qq pair 

with the definitions 

7 = , - £ —. (1.2) 

This form is not very sensitive to the orientation of the axis of the pair as compared with 
its velocity. In particular, the average value of the maximal length over the corresponding 
angle is, under reasonable assumptions on the angle distribution of the produced qq pairs, 
basically of the same form as well. In the above, N is the number of colors, T is the 
plasma temperature and J is the R-charge density. Also, ^ is the dimensionless R- 
charge parameter (the overall numerical value in its definition will be fixed later), while 
A = A(,^) is a dimensionless function of £ ranging between zero and an upper bound 
of order one. Finally, q, with i = 1,2, 3, are numerical parameters which, although in 
our models they have precisely determined values, relatively small compared to unity, in 
more realistic experimental situations they should be fitted to data. The leading term in 
Eq. (1.1), which dominates in the ultrarelativistic (large-7) limit, was already suggested 
by the computation of [21]. However, we shall see that in our models it makes sense to 
include the other two terms since they become important for not too extremely relativistic 
velocities, i.e. for v < 0.8. The function T{v,^) has a very weak functional dependence 
and can be practically taken to be constant. The function G{^) has a behavior of the 
form G{^) = 1 + O^^"^) for small ^-values, but can otherwise have cither weak or strong 
dependence on the parameter ^, depending on the model. In particular, in cases with 
^ ^ 1, we have the behavior G{^) ~ signifying that the screening length becomes 

temperature-independent and its scale is set by the R-charge density J. In addition, a 
dependence proportional to N'^^^ arises. For experimentally-motivated studies however, 
one may compare ratios of maximal screening lengths at different velocities so that any 
strong dependence on the chemical potential cancels out. We note that in [24] it has been 
shown that the leading velocity dependence of the screening length is of the form {l—v'^Y, 
where the exponent = ^ in A/" = 4 theories, but is lower than that in certain cases with 
less super symmetry. Based on that we may further generalize (1.1) and replace in the 
exponent of the leading order term the 1/2 by the parameter 2v. We do not believe that 
this generalization will affect the form of the correction terms, as they appear essentially 
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from the analyticity properties of integrals below. The same replacement can be done 
for the analogous formulae below. 

In the rest frame of the plasma, we should also take into account the Lorentz contraction 
of the component of the length along the direction of the motion. In this case the result 
is sensitive to the actual distribution probability for the angle between the qq-pair axis 
and its velocity vector. The most reasonable assumption is that, that this probability 
is uniform in the plasma rest frame. Then, based on our computations, we suggest the 
following law 

Plasma rest frame(l) : L^ax — ^{v,^) ^ ( 1 + di— ' ^"^ ~'~ ^^"^ 



(1.3) 

where again we have averaged over the angle formed by the gg-pair's axis and its velocity 
vector. The numerical parameters a and di, i = 1,2,3 are to be fitted to possible 
experimental results. Note that, compared to (1.1), we have a suppression factor I/7 as 
as well as the presence of a In 7 term in the overall coefficient. As we will show, these 
are due to the dominance, in taking the average, of pairs with their axis almost parallel 
to their velocity vector. 

A second possibility, less likely in our opinion, is that the angle distribution be uniform 
in the pair's rest frame. In this case, instead of (1.3) we have 

Plasma rest frame(2) : I/max = ^{v, i) —7^ —77^ ( I + di — + ~'~ + ^4 In 7 



(1.4) 

where again we have averaged over the angle. Note that in this case we have no extra 
suppression in the length compared to (1.1), with the exception of the appearance of In 7 
in a subleading term. Hence in this case we have a dominance of angles near © = 7r/2, 
at which the length suffers no Lorentz contraction. 

The organization of this paper is as follows. In section 2 we set up the general formalism 
for computing the interaction potential of moving qq pairs in general backgrounds. In 
section 3 we review the zero R-charge case and in section 4 we compute the maximal 
screening length for more general situations having non-zero R-charge using rotating D3- 
brane solutions. Finally, in section 5 we consider issues concerning the screening length 
in the plasma rest frame versus the quark pair rest frame as well as their average values. 
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2 Moving qq pairs in the Af = 4 plasma 



In the framework of AdS/CFT, Wilson-loop calculations for a qq pair moving with ve- 
locity V relative to the plasma can be carried our either in the rest frame of the plasma 
where one considers a string whose endpoints move with velocity v or, equivalently, in the 
rest frame of the qq pair, where one considers a string with fixed endpoints in a boosted 
background. In what follows, we will employ the second approach. 

To describe the method, we begin with a general ten-dimensional metric of the form 

ds^ = Gttdf + G^^dx^ + Gyydy^ + Grrdr^ H , (2.1) 

where the metric components do not depend on t, x and y and where the ellipses denote 
the remaining terms involving all other variables. These could be either cyclic variables 
similar to x and y or angular variables that can be set to some constant values consistent 
with their equations of motion. Orienting the axes so that the qq pair moves along the 
negative x direction, and applying a boost in this direction, we obtain the metric 

ds"^ = Gttdt'^ + G^cxdx^ + 2Gtxdtdx + Gyydy"^ + Grrdr'^ H , (2.2) 

where 

Gtt = I'iGu + v'^G,,) , G,, = ^'{v'^Gu + G,,) , G*, = -i''v{Gu + G..) ■ (2.3) 

The separation length and energy of a qq pair of arbitrary orientation with respect to 
the plasma velocity is found by considering a Wilson loop on a rectangular contour C 
with one side of length T along the t direction and one side of length L along a spatial 
direction, taken for definiteness to be parallel to the xy plane and tilted at an angle 
relative to the x axis, as shown in Fig. 1. Then, according to the gauge/gravity 
correspondence, the interaction potential energy E of the qq pair, as measured by the 
Wilson loop expectation {W{C)), is given by [2] 

e-^-^^ = {W{C)) = e'^t^l , (2.4) 

where 

S[C] = j dTda^J-detg^f) , g^p = G^^daX^dpx" , (2.5) 

is the Nambu-Goto action for a string propagating in the boosted supergravity back- 
ground, whose endpoints trace the contour C. We note that we will use Minkowski signa- 
ture throughout the paper. Fixing reparametrization invariance by taking (r, a) = {t, y), 
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Figure 1: Wilson loop for the energy computation of a qq pair in a moving plasma (shaded 
rectangle). 



assuming translational invariance along the t direction and setting all angles to constants 
(to specific values consistent with the equations of motion), the embedding of the string 
is described by the functions 



r = u{y) , X = x{y) 



rest = const. 



and the boundary condition for the string endpoints becomes^ 



u ( ±— sin 6 
2 



oo 



X ±— sin B 
2 



±— cos 9 . 
2 



The Nambu-Goto action is then written in the form 

S[C] = -^J dy^fy{u) + Uu)x'^ + g{u) 
where the prime denotes a derivative with respect to y and^ 



/2 



-^^Grr{Gu + V^G, 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



^Strictly speaking, the above embedding is not valid at 8 = 0, in which case one has to interchange 
the roles of x and y i.e. take (r, a) = it, x) and set r = u{x) and y = y{x). However, the final formulas 
for that case, as quoted after Eq. (2.14), can be obtained as limiting cases with the present embedding. 

■^Note in passing that, in the limit t; ^ 1 (7 00) we get the action used in [14] in supergravity 
computations of the jet quenching parameter. In taking this limit the extra 7 factor is absorbed in the 
time dilation of the long side of the Wilson loop. 
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Since actions of the precise form (2.8) have been encountered in the static Wilson-loop 
calculations for a;' = in [4] (see section 1) and [25] (see section 5, where x played the 
role of an angular variable), we briefly review the results without entering into the inter- 
mediate computational details. Independence of the Lagrangian from y and x leads to 
two first-order equations expressing the conservation of the "energy" and the momentum 
TTx conjugate to x{y). Integrating these equations, one obtains [25] 

du 



Lsine = 2yr^/V2(Ko) / y- 

Juo jy 




Lease = 27r./V^(«o) r^^/^, (2.10) 

Juo Jx{u) y F{U) 

where the integration constant Uq is the value where u{y) develops a minimum. Substi- 
tuting them in the Nambu-Goto action, one finds the qq potential energy 




E=zl du,lf^^-E,. (2.11) 

where Eq is the subtraction term corresponding to the energy of two disconnected world- 
sheets, to be discussed below. In the above, the function F{u) is defined by^ 

and satisfies F{uq) — 0. The system of Eqs. (2.10) implicitly determines the parameters 
Uq and TTj; in terms of the separation length L and the angle ©, the latter clearly only via 
the combinations cos^ © and sin^ ©. Substituting the resulting values into (2.11) gives 
the potential energy for the same parameters. Note that Eqs. (2.10)-(2.12) are invariant 
under the symmetry operation © ^ | — © followed by the exchange of the x and y 
indices and the renaming tt^ ^ 1 — tt^. 

To avoid unnecessary complications, we will restrict to the cases where the spatial side 
of the Wilson loop is either perpendicular (© = |) or parallel (© = 0) to the plasma 
velocity and only discuss the general case in section 5. In the first case, we have tTj; = 
and so the length and energy read [4] 



= Vl'\u,) rduJ , (2.13) 



and 




duj fff'' ^ -E„. (2.14) 



'Compared to [25], this function has been rescaled by a factor of fy{u). 
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while in the second case we have iix = 1 and the length and energy are given by the 
above formulas with the replacement fy ^ fx- 

At this point, some clarifying comments are in order. A first issue is how much the 
string may stretch into the radial direction, i.e. which is the minimum allowed value 
for the turning point Uq. To examine this issue, we note that in the geometries under 
consideration there appear several special values of the radius, namely (i) the horizon 
radius uh where g{u) blows up, (ii) the radius Ug where fx{u) becomes zero, which in our 
examples, depending on the angle, coincides with the maximum or minimum values of the 
ergosphere (which, in the absence of rotation, coincides with Uh) and (iii) the velocity- 
dependent radius where g{u) and fy{u) become zero^ (which, in the zero- velocity limit, 
coincides with uh)- In our examples, these radii always satisfy uh ^ Ug ^ Uy. While it 
is quite obvious that taking uq > Uy is always consistent, some care is needed in order to 
examine whether we can take uo < u^. To this end, we distinguish between two cases. 
When © 7^ 0, the first of Eqs. (2.10) is always relevant and the factor y^Tjuo) in that 
equation becomes imaginary for uq < u^, therefore excluding that possibility. On the 
other hand, when © = 0, the first of Eqs. (2.10) becomes irrelevant and the behavior of 
the string is solely determined by the function x(u) with x'(u) ~ \/ fxiuo) — ^^^^ — If 
we take uq < then the string will first reach the radius where x'{u) becomes zero so 
that the string develops a cusp and cannot be further extended towards uq. If one wishes 
to include such configurations, one must thus cut off the region uq < u^hj replacing 
the lower limit of the integrations from uq to Uy for all Uq < u^. However, since (i) string 
configurations with cusps raise questions about the validity of the Nambu-Goto approach 
since the string curvature diverges there, (ii) the maximal value of L{uq), which is the 
main quantity of interest in this paper, corresponds to a value Uq of uq that is always 
larger than u^, and (iii) the region uq < Uq where cusps may appear corresponds to a 
branch of the solution (with the same prescribed boundary conditions (2.7)) which is 
energetically unfavorable and thus presumably unstable, we will altogether ignore the 
presence of cusps as they are irrelevant for our screening problem. 

Another issue refers to the choice of subtraction term Eq in (2.11). In the zero-velocity 
case, the subtraction term is unique and represents two straight strings stretching from 

^The significance of the radius Uy has been pointed out in [8] for one-quark configurations and in 
[21, 22, 23] for two-quark configurations. Roughly, u = Uy acts as a hypersurface for the boosted metric, 
beyond which the string cannot penetrate unless a special choice of integration constant is employed in 
the one-quark case [8]. 
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infinity to the horizon. In the case of nonzero velocity, however, the above configuration 
is no longer a solution to the equations of motionand one must find an appropriate 
subtraction term that (i) removes the divergence, (ii) is consistent with the equations of 
motion, (iii) is independent of the separation length and (iv) reduces in the zero- velocity 
limit to two straight strings stretching to the horizon. These criteria are satisfied by any 
configuration of two disconnected strings interpolating between the drag configuration 
of [8] (and its generalizations [9] for the case of rotating branes), which represents a 
curved string stretching up to the horizon uh and a configuration of one straight string 
stretching up to the radius Uy ^ uh- Note that in order for the strings to end at a 
radius larger than uh, one needs to assume the presence of a second flavor brane at that 
radius, in which case each string should be regarded as dual to a meson composed out 
of a heavy and a light quark (like a bottom and an up or a down quark). Be that as 
it may, there seems to exist an infinite family of subtraction terms [26] consistent with 
the above requirements. The maximum subtraction energy [22] corresponds to the drag 
configuration and is given by 

1 r°° 

Eo = - du , (2.15) 

JuH 

while the minimum subtraction energy corresponds to the straight-string configuration 
and reads 

Eo^- du^f^) . (2.16) 

Among the possible subtraction terms, the one corresponding to the drag configuration 
seems to be the most sensible one, as it is the only solution available in the absence of 
an extra flavor brane. 



3 Screening in the zero— R-charge case 



In order to establish our notation and present some extra details, we first review the zero- 
R-charge case, examined in [21] and also in [22]. In this case, the background geometry 
has metric 



l-^i^\df + dx' + dy' + dz' 



2^2 



+ 



dr^ + R^dnl , 



where the horizon is located at r = /x and the Hawking temperature is T = 
we have for the various functions 



fy{u) = -^{u - 7 ) , fx{u) = -^{u -n) , g(u) = 



(3.1) 
Then 

(3.2) 
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3.1 Motion perpendicular to the qq-pair axis 



Let us first consider the © = | case corresponding to plasma velocities perpendicular 
to the axis of the qq pair. Substituting fy{u) and g{u) from (3.2) into Eqs. (2.13) and 
(2.14), introducing the dimensionless length and energy parameters 

1 



e = ttTl , 



£ — 



E 



(3.3) 



and trading u and uq for the dimensionless variables y — ^ and yo — ^^ s/j, we obtain 
for the length 



[v,yo) 



yt 



dfj 



2/0 



Viy'-^W-yl) 



yl 



2-^1 I n> 7' 7' 77 



(3.4) 



r(l/4)2 y-^^ V2'4'4'|/o^ 

where 2Fi{a,b,c,x) denotes the standard hypergeometric function. For the energy, use 
of the maximum subtraction (2.15) leads to the result 



^{'v,yo) 



PCX 

Jyo 



dy 



y -7 



^{y^-l){y^- 

yo 2F1 



- 1 



1111 

4'2'4'^ 



r(l/4)2 

while use of the minimum subtraction (2.16) yields 



- {yo - 1) 



yt 





5 




V2'4' 


4' 


yt). 



+ 1 , (3.5) 



£(^^,yo) 



dy 



yo 



r 



^{y'-m'-yt) 



r(i/4)2 



1111 

1/0 I "i' 2' 4'^ 




(3.6) 



-2x/7 2F1 



1 5 1 



2' 4'4'72_ 

Note that, according to the discussion of section 2, no cusps occur and that the two 
expressions above differ only by the last term which is velocity dependent, but totally in- 
dependent of the parameter uq. As a result the two expressions differ by an ^-independent 
term and the force between the quarks, given by minus the derivative of the energy with 

respect to the separation length, is the same no matter which expression we use. In the 
static limit in which 7=1, using properties of the hypergeometric functions, the expres- 
sions (3.5) and (3.6) are found to coincide. Also for £ — > 0, both go to the conformal 
result of [2]. 
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(a) (b) (c) 

Figure 2: (a,b) Potential energy e plotted as a function of the separation £ for the 
values f = 0, 0.5, 0.7, 0.9 and 0.99 (right to left) using the maximum and minimum 
subtraction prescriptions respectively, (c) Plots of various approximations to dmaxiv), 
namely the result of numerical computation (solid), the approximation (1.1) (dashed) 
and its leading-order form (dotted). 

Fig. 2(a) shows the dimensionless energy e as a function of the dimensionless separation 
length ^ for various values of the velocity v. We see that e is a double-valued function 
of £ whose upper branch corresponds to a presumably metastable, at best, configuration 
and is to be discarded. As in the static case with f = [3], there is a maximal length 
^max{v) beyond which there is complete screening of the potential since the only allowed 
configuration is the zero-energy configuration of two disconnected worldsheets. Its nu- 
merical value for f = is i'max(O) = 0.869 was given in [4, 21]. For ultrahigh velocities 
it goes to zero as (-raa^iv ^ 1) = 0.7437^^/^ as first noted in [21]. There are various 
ways to see that. In particular, assuming that Uq ^ 1, from (3.4) and the fact that the 
hypergeometric function becomes unity at the origin, we see that £ ~ '^■^'ilil ^^"^"'^ . 
Maximization of this with respect to yo gives ?/g = v^7, which upon substitution gives 



'-max 



f ^ 1) = 33/4r((/4-)2 7 ^^'^ — 0.7437 as stated above. 



Besides the maximal length, there also exists a value isiv) of the length [21, 22] (its value 
for f = is 4(0) = 0.754 [4]) after which the lower branch of e takes positive values. 
The region beyond this point is also to be discarded, since the free configuration becomes 
energetically favored there, and hence this length could be also used as a definition of 
the screening length. Using the regularization (3.5) for the energy there is a velocity 
Vc — 0.447 for which the maximum length and is coincide [22] and after which there is 
no distinction between them. Then, as shown in Fig. 2(a), for small velocities the length 
is{v) differs from the maximal one by about ten percent and, as the velocity is increased. 
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the difference becomes smaller and vanishes for the above value of the velocity Vc- Since 
£s practically coincides with ^raax for v > 0.35, and since its study necessitates examining 
both e and ^ thus requiring more involved numerics, in this paper we will define the 
screening length as the maximal length allowed. 

The above results for £max('i') served as a motivation for the authors of [21] to suggest 
that the velocity dependence of the maximal length is given by the empirical formula 

ir.Uv)-Hv){l-v^Y/\ (3.7) 

where J-{v) is a function with relatively weak dependence on v, with .^(0) = 0.869 and 
!F{1) — 0.743. With an analysis that we will perform in a more general setting in the 
next section, we will find an improved version of (3.7) 

Cax(^) = - [1 + C2(l - V^f'^] , (3.8) 

which is of the form (1.1). For our model the numerical coefficient is computed in more 
general situations below and takes the value C2 = 0.100. In that case T{v) is indeed a 
practically flat function of v since J^{Q) = 0.790 and J^(l) = 0.743. Plots of imaxiv) in 
the various approximations considered are shown in Fig. 2(c). 



3.2 Motion parallel to the gg-pair axis 

Let us consider next the © = case, corresponding to plasma velocities along the axis of 
the qq pair. Substituting fx{u) and g{u) in the corresponding equations and employing 
the same rescalings as before, we find for the length 



£{v,yo) = 2^1/4-1 [ 



dy ly^- 7^ 



,0 y'-i\y'-yl ' 



r(l/4)2 y3 i^4' 2' '4'yo^'yo^ 

where Fi{a,bi,b2,c,x,y) is the Appell hypergeometric function. For the energy, use of 
the maximum subtraction (2.15) and the minimum subtraction (2.16) leads to 




£(v,yo) = / dy\A^ j-1 -(yo-1), 



r(l/4)2^"" 1 1^ 2' 4' 4' 



.,2 



yo2Fi{--,--,-,^]+l. (3.10) 
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and 



r(i/4)2 



r 



yo 2F1 { 




(3.11) 



respectively. Note that, according to the discussion of section 2, one may take yo < 
Vv = \/7 in which case there may occur cusps. However, the branch where cusps occur 
is irrelevant for our problem. This subtlety aside, the qualitative behavior discussed in 
the previous subsection remains virtually unchanged, with only small deviations from 
the © = f case. For instance, the numerical coefficient for the i> — > 1 behavior changes 
to 0.837 but nevertheless the velocity dependence remains of the same form. For this 
reason, in the rest of this paper we will consider in detail only the © = f case. 



4 Screening in the case of nonzero R-charges 

In light of the above results, a natural question arising is whether the (1 — v'^Y^'^ velocity 
scaling law is universal in more generalized settings, one of which is the Af — A theory with 
non- vanishing R-charges. In what follows, we will briefly review the relevant supergravity 
backgrounds and then we will address the questions (i) whether this velocity dependence 
persists in the presence of R-charges and (ii) what is the nature of the dependence of the 
maximal length on the R-charge density. 



4.1 Non-extremal rotating branes 

According to AdS/CFT, the gravity duals to A/" = 4, SU{N) SYM theory at finite 
temperature and R-charge are given by non-extremal rotating D3-brane solutions found 
in [27, 28]. In the conventions of [29] which we here follow, the field-theory limit of these 
solutions is characterized by the non-extremality parameter /i and the angular momentum 
parameters Oj, i — 1, 2, 3, which are related to three R-charges (or chemical potentials) 
in the gauge theory. Below we review two special cases of these solutions, corresponding 
to two equal nonzero angular momenta and one nonzero angular momentum, for which 
we state the relations between the supergravity and gauge-theory parameters as well as 
the restrictions on the physical ranges of these parameters as dictated by thermodynamic 
stability [30]. For more details on the calculations, the reader is referred to [14]. 
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4.1.1 Two equal nonzero angular momenta 



We first examine the case of two equal nonzero angular momenta, which we can take as 
ci2 = 03 = ''^o- The metric is given by 



dr + dx^ + dy^ + dz^ 



(r4-/x4)(r2-r2) 



(r^ - rl cos" e)de'' + cos" ^dQ^ + (r" - rl) sin^ ^dc/,^ 

2 

cos^ 6'(sin2 + cos^ ^'#3)! , 



(4.1) 



where 



and d^l — dt/j^ + sin^ ipd(l)l + cos^ ipdcpl is the -S"^ hue element. The horizon is located at 



(4.3) 



On the gravity side, the above solution is characterized by the two parameters and tq. 
On the dual gauge-theory side, it is most natural to use the Hawking temperature T and 
the angular momentum density J, which correspond to the temperature and R-charge 
density of the gauge theory i.e. to the canonical-ensemble thermodynamic variables. On 
both sides, it is convenient to trade tq and J for the dimensionless parameters 

6^37 



A 



The two sets of parameters (//, A) and (T, ^) are related by 

Q 



(4.4) 



(4.5) 



where F{^) is the function 



Q{i) 



{i + ^^TeY'' - (e + yiTer"^ 

a/3 



(4.6) 



Substituting this into (4.5) wc obtain the explicit expression for A = A(.^). The constraints 
imposed by thermodynamic stability on the parameters A, ^ and Q as follows 



O^A^l, 0^C<oo, 0^Q<oo. 



(4.7) 
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4.1.2 One nonzero angular momentum 



We next examine the case of only one nonzero angular momentum, which we can take 
as ai = Tq. The metric is given by 



- 1- 



+ dx^ + dy^ + dz^ 



+ 



2^2 



1/2 



(r^ + cos^ ^)cie2 + cos^ ^dQ^ + (r^ + r^) sin^ ^d^^ 



(4.8) 



2^sin2 



where 



j.2(j.2 _|_ ^2 j,Qg2 Q-\^ 

and the line element is as before. The horizon is located at r = r^^ with 



(4.9) 



(4.10) 



As before, the solution is characterized by /i and vq on the gravity side and by T and J 
on the dual gauge-theory side. Introducing again dimensionless variables according to 

2V2J 



(4.11) 



we find that the two sets of parameters are related by 



1/4 



A 



V2{Q - 1)1/2 



(4.12) 



Ql/4(2-Q)l/4 ' 

where Q{^) is the lower branch of the solution to the equation Q'^{Q — l){Q — 2)+^'^ — 0, 
which admits a perturbative expansion, near ^ = 0, as Q — 1 + 0{^'^). Again, one 
may substitute this into (4.12) to obtain a perturbative expansion for A = A(^). Finally, 
thermodynamic stability constrains A, ^ and Q as follows 



0^A^Ao~1.46, O^^^emax^l.33, 1 ^ g ^ Q^a. ^ 1-73 . (4.13) 



4.2 Velocity and R-charge dependence of the maximal length 

After the above preliminaries, we are ready to examine the problem of screening in the 
presence of nonzero R-charges. The new parameter entering the problem in this case is the 
supergravity angular momentum parameter A or, equivalently, the gauge-theory R-charge 
density ^. Therefore, the dimensionless length and energy have the form i{v,^,yo) and 
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e{v,C,,yo), while the maximal length is of the form ima.x{v,C)- In the rest of this section 
we focus on the maximal length and on its dependence on the velocity and R-charge. 

As a general remark we note that the leading high-velocity behavior of the screening 
length in the presence of R-charges will be the same as that in their absence, since it 
occurs for large values of the parameter uq proportional to 7^/^, for which the geometry 
behaves as if it was AdS^ x S^. This was also noted in [24] from a five-dimensional 
perspective^ and is apparent when the comparison is done at constant energy density 
e = above extremallity, in which case the (dimensionful) maximal length is of 

the form Lmax ~ ^^^^'J ^ ~ "ttt"! with the same coefficient as in the zero R-charge 
case. However, in this paper we compare quantities at constant temperature so that the 
proportionality factor is apparently different and depends on the R-charges. 

4.2.1 Two equal nonzero angular momenta 

For this case, the functions fy{u) and g{u) are given by 



where for consistency with the corresponding equations of motion the angular variable 6 
must be either or tt/2. The expressions for the dimensionless length £ are as follows. 

• For the case with ^ = 0, after a change of integration variable, the dimensionless 
separation length takes the form^ 



At this point, we are ready to justify for this model the proposed empirical formula 
(1.1) for the screening length. Since maximization of (4.15) with respect to yo cannot be 

^We emphasize that results from the Nambu-Goto actions with a ten-dimensional background metric 
are generically different than those obtained with the background metric replaced by the five- dimensional 

one corresponding to a lower dimensional supergravity theory. A recent such example is the computation 
of the jet quenching parameter in the presence of R-charges. The results of [14, 15], where a ten- 
dimensional approach was followed agree only qualitatively with those of [16], whose approach is a 
five-dimensional one. 

^In the limiting case A = 1 (corresponding to ^ ^ 00), the expression (4.15) and the similar one 
(4.18) below can be evaluated exactly in terms of complete elliptic integrals. However, no such exact 
formulas are available for the integrals (4.23) and (4.26) further below. 




^v, ^, yo) 



Vyp - A^i/g - 7^ 




(4.15) 
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(a) (b) (c) 

Figure 3: (a), (b) Plots of imaxiv) for the case of two angular momenta at 6* = | (with 
similar plots for 6 = 0) with A = 0.5 and A = 0.95 respectively; The solid, dashed and 
dotted lines are as in Fig. 2(c). (c) Plot of the function G{^) for the case of two angular 
momenta. 

carried out analytically, we employ a perturbative method that proceeds as follows. We 
first assume that the maximum, for large enough 7, is at some value of yo of leading order 
7^/^, to be justified by the result, and then we expand the integrand into inverse powers of 
Uq up to second order. The resulting expression can then be analytically maximized with 
respect to yo and then the result for £max follows by substitution. We find an expression 
of the form (1.1) with 

ci = -0.10525 , C2 = 0.100 , C3 = -0.01507 , (4.16) 

where we preferred to give, to a certain accuracy, the numerical values of the various 
coefficients, instead of presenting their explicit expressions in terms of gamma functions. 
For jF(t>,^) we give its numerical values in the four corners in the (f,^)-plane 

:F(0,0) = 1.06 , J^(0, 00) = 0.96 , J^(1,0) = 1, J^(l,cx)) = l, (4.17) 

showing that it is indeed a very slowly varying function. 

• For the case with 6 = 7r/2, we have the expression 

i{v,^,yo) = G{0 r ^ =. (4.18) 

Performing a similar analysis as before we find that the various coefficients entering (1.1) 
are given by 

ci = 0.2106 , C2 = 0.100 , C3 = 0.0898 . (4.19) 
For J-'{v,^) we give its numerical values in the four corners in the (f,.^)-plane 

J^(0,0) = 1.06 , J^(0,oo) = 2.11 , J^(1,0) = 1, J^(l,oo) = l, (4.20) 
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showing that it is a slowly-varying function except for the case of low velocity and ex- 
tremely large R-charge. However, the region of the (f , ^) parameter space where these 
deviations from unity appear is very small, as can be easily verified by an explicit plot 
of^KO- 

In both cases the function G(^) is given by 

G'(e) = (l + F2)-V2 = ^r3A2 . (4.21) 

It is plotted in Fig. 3(c) and it has a strong dependence on the dimensionless R-charge 

density ^. 

4.2.2 One nonzero angular momentum 

Now, the functions fy{u) and g{u) are given by 

and again the angular variable 6 — or 6 = ^ for consistency. The length integrals are 
given below. 

• For the case with ^ = after a change of variables we find that 



dz 



An analysis similar to before gives for the numerical factors in (1.1) that 



(4.23) 



ci = 0.10525 , C2 = 0.100 , C3 = -0.01507 . (4.24) 

The numerical values of T{y^ ^) in the four corners of the {y^ ^)-plane are 

JP(0, 0) = 1.06 , JP(0, 1.33) = 1.05 , 0) = 1 , T{\, 1.33) = 1 , (4.25) 

showing again that it is indeed a very slowly varying function. 

• For the case with ^ = 7r/2 we obtain after a similar variable change that 

Now the numerical factor in (1.1) are 

ci = -0.2106 , C2 = 0.100 , C3 = 0.0898 , (4.27) 
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(a) (b) (c) 

Figure 4: (a), (b) Plots of £max(^) for the case of one angular momentum at ^ = (with 
similar plots for 6* = |) with A = 0.7 and A = 1.4 respectively; the solid, dashed and 
dotted lines are as in Fig. 2(c). (c) Plot of the function G{^) for the case of one angular 
momentum. 

while the numerical values of J^{v, ^) in the four corners of the {v, ^)-plane are 

J^(0,0) = 1.06 , J^(0, 1.33) = 0.74 , J'(1,0) = 1, J^(l, 1.33) = 1 , (4.28) 

showing that it is a slowly- varying function except for the case of low velocity with large 
R-charge when further corrections become important. 

Now we have for the function G{^) the expression 

I / I \ 

G(0 = F-^/^(2 - F)-/^ = Vi + T(r + T-Y) • ^^-29) 

This is depicted in Fig. 4(c) and shows a very weak dependence on the dimensionless 
R-charge density ^. 

One comment is in order here regarding the case of parallel motion having = 0. In 
that case the form of the maximal screening length is the same with slightly different 
coefficients and function JF. The function G{C,) is the same and is independent of the 
angle 0, a fact that facilitates the computation of average values of the maximal screening 
lengths in the next section. 

5 Plasma rest frame and averages 

So far we have not considered the most general situation where the gg-pair axis makes an 
arbitrary angle G with the velocity vector. A detailed investigation of this problem can be 
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carried out cither by numerical analysis of the full equations (2.10) or by an approximate 
analytic method that we will present here. We have already noted that the maximum 
separation length is a function of through the combinations cos^ and sin^ 0. Since 
its values for the extreme case = and = 7r/2, do not differ very much, it is a good 
approximation for the screening length in the rest frame of the pair to write 

-f'max,qq(0) = siu^ + L\\ COS^ , (5.1) 

as it reduces to the above mentioned extreme cases, with the maximum screening lengths 
denoted by L± and L||, in an obvious notation. It is of interest to know the average 
vahie of this length over all angles. To do that we have to know the probability for a qq 
pair to be moving at an angle, formed by its axis and its velocity vector, between and 
+ dQ. The most reasonable assumption is that this probability is uniform, but in the 
plasma rest frame. Denoting in the plasma rest frame by 0' the angle corresponding to 
the angle in the qq rest frame, we have, due to length contraction in the direction x 
parallel to the velocity vector, the relation 

tan0' = 7tan0, (5.2) 

that is the length is not only contracted but also rotated. From these the useful relations 

9^/ 7"^cos^0 . sin^0 . 

7 ^ cos^ + sm 7 ^ COS"' + sm 

and their inverses 

2n— ",-eos-0' • 2(-v_ sin-0' 

" 72cos2 0' + sin2 0' ' '''' ^ " 72cos2 0' + sin^0' ' ^ ' 

follow. Then the probability distribution is 

Ppla.ma(e')d0' = -dQ' = —— „^ , . = P,,-(0)d0 . (5.5) 

TT 7r7 7-2 cos^ + sm'' 

Then the average value of the maximal screening length in the gg-pair rest frame is 

Lmpk^.an / C?0 

Jo 



max,qq I o o , ■ 1 /-\ 

7r7 Jo 1 cos^ + sm^ 

L„ + (L.-L„)(i-iUo(A) . (5.6) 



1 + 7 ~" ' Vt 7^/ V7^ 

This is essentially of the form (1.1)*^ and for large enough 7 it is basically controlled by 
the maximal screening length L\\ along the pair's velocity vector. The reason for this 



^ There are two additional terms appearing in the parenthesis of the form I/7 and 
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is the fact that the probabiUty (5.5), although uniform in the plasma rest frame, it is 
peaked at = with a width of order I/7, in the gg-pair rest frame. 

In the plasma rest frame the expression corresponding to (5.1) differs from the above 
in that we have to take into account the Lorentz contraction of the component of the 
separation length vector parallel to the velocity. We find 



^m£«,piasma(©') = {L± siu^ 9 + L|| cos^ 6) -^7-2 cos^ 6 + siu^ 6 
L^sin2e' + L||72cos2e' 



with average 



(72cos2e' + sin2eO=^/2 ' 



-^max, plasma / dQ 



(5.7) 



p7] 

Jo 

[L^ + L|| ln(4e-S)] + O (^1^ . (5.8) 



^1 Jo ^7-2 cos2 e + sin^ e 

_ 2 
777 

We see that there is a I/7 suppression of the maximal screening length, the reason being, 
as explained before, that in the pair's rest frame the production is peaked near = 0. 
The presence of the In 7 term in the overall coefficient is due to the contribution of the 
integral in the region near = 0. The last expression in the basis of our proposal for a 
maximum screening length in the plasma rest frame in (1.3). 

We note for completeness that, if we take the angle distribution to be uniform in the 
g^-pair rest frame, then the average value of the maximum length in the rest frame of 
the pair would be, instead of (5.6) 

- _ L± + L|| 

-'^max,qq — ^ ' \'-'-'^J 

again of the form (1.1). Also the corresponding average in the plasma rest frame would 
be 

2 f^^^ I 

-^max,plasma = - / dQ {L± siu^ + L|| COS^ 0) J 7"^ COS^ + siu^ 

Jo 



1 

2 
2 

37r 



,Fi(-l,l,2,^^) + L|| ,Fi(-l,^,2,^2) 
(2L^ + L|,) + J-^ (2L^ - 5L|| + GL,, ln(47)) + O (^-^^ 



(5.10) 



The last expression is the basis of our second possibility, less likely in our opinion, for a 
maximum screening length in the plasma rest frame in (1.4). 
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